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The transcriptional regulatory protein Fnr, acts as an intracellular redox sensor regulating a wide range of 
genes in response to changes in oxygen levels. Genome sequencing of Herbaspirillum seropedicae SmRl 
revealed the presence of three fnr-like genes. In this study we have constructed single, double and triple fnr 
deletion mutant strains of H. seropedicae. Transcriptional profiling in combination with expression data 
from reporter fusions, together with spectroscopic analysis, demonstrates that the Fnrl and Fnr3 proteins 
not only regulate expression of the cbb 3 -type respiratory oxidase, but also control the cytochrome content 
and other component complexes required for the cytochrome c-based electron transport pathway. 
Accordingly, in the absence of the three Fnr paralogs, growth is restricted at low oxygen tensions and 
nitrogenase activity is impaired. Our results suggest that the H. seropedicae Fnr proteins are major players 
in regulating the composition of the electron transport chain in response to prevailing oxygen 
concentrations. 



■w- -w- . seropedicae is an endophytic diazotroph belonging to the Betaproteobacteria that can fix nitrogen 
m m under micro-oxic and nitrogen limiting conditions. It is found in association with different crops such as 
■M. J. rice, maize, sugar cane and sorghum 13 . H. seropedicae SmRl is an aerobic bacterium which has a 
branched respiratory chain comprising different types of terminal oxidases, which potentially could allow the 
bacteria to exploit respiratory flexibility and survive under microaerobic conditions 4 . Genome sequencing 
revealed genes coding for three Fnr-like proteins in H. seropedicae SmRl 4 . In many organisms the Fnr protein, 
which belongs to the CRP-FNR family of transcriptional regulators 5 , acts as a positive or negative regulator of 
genes required for the metabolic switch in response to 0 2 levels 6 . The Fnr protein of Escherichia coli contains an 
N-terminal sensory domain, which binds an oxygen-labile [4Fe-4S] 2+ cluster under oxygen limiting conditions 7 
and a C-terminal DNA-binding domain, which recognizes a partially palindromic sequence called the Fnr-box or 
anaerobox, TTGAT-N 4 -ATCAA 5,S . 

Various Fnr-related transcriptional regulators of the CRP-FNR family have been reported to be involved in 
biological nitrogen fixation. The Fnr protein from Klebsiella pneumoniae is required to relieve inhibition of NifA 
activity by its partner regulatory protein NifL under anaerobic conditions 9 . In Rhizobium leguminosarium 
UPM791 FnrN is responsible for the expression of the high affinity oxidase encoded by fixNOQP that supports 
growth under microaerobic conditions and is essential for nitrogen fixation 10 . Similarly the FixK 2 protein is also 
essential for nitrogen fixation in Bradyrhizobium japonicum and Sinorhizobium meliloti 11 ' 13 '. Representatives of 
the CRP-FNR family are also known to act negatively in repressing genes related to nitrogen fixation, such as the 
FixKi protein from B. japonicum 12 . 

As the H. seropedicae genome encodes three Fnr-like proteins, we were interested to determine the potential 
involvement of these three Fnr homologs in nitrogen fixation and in the control of gene expression in response to 
oxygen limitation. Several representatives of the Betaproteobacteria encode more than one Fnr-like protein in 
their genome. For example, Burkholderia pseudomalei 1710b and Herminiimonas arsenicoxydans have genes 
coding for two Fnr-like proteins, whereas Cupriavidus metalidurans CH34 and Ralstonia eutropha H16 encode 
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three and five Fnr-like proteins, respectively. However, to date the 
functions of these Fnr-like paralogs have not been determined. 

In the current study, we sought to attribute function to the three 
Fnr-like proteins found in H. seropedicae and in particular to exam- 
ine their role in the regulation of electron transport chain composi- 
tion. We demonstrate that deletion of all three fnr alleles results in a 
growth phenotype under microaerobic conditions, implying that Fnr 
proteins may be involved in controlling the expression of respiratory 
oxidases. By comparing the transcription profiles of the wild-type 
and triple fnr mutant strains and performing further gene expression 
and biochemical analyses, we observe that the Fnr proteins not only 
activate genes required for expression and activity of the high affinity 
cbb 3 -type oxidase, but also have a major influence on the regulation 
of the cytochrome bc x complex and on cytochrome c biogenesis. This 
suggests that the H. seropedicae Fnr proteins facilitate distribution of 
the electron flux through cytochrome carriers and the heme-copper 
oxidase branch of the respiratory chain to increase coupling effi- 
ciency under oxygen -limiting conditions. 

Results 

H. seropedicae encodes three proteins that share homology with 
E.coli Fnr. The H. seropedicae SmRl genome 4 contains three genes 



encoding homologs of the Fnr protein, which we designate as fnrl, 
(Locus Tag: Hsero_3197; Ref seq: YP_003776587.1), fnr2 (Locus 
Tag: Hsero_2381; Ref Seq: YP_003775788.1) and fnr3 (Locus Tag: 
Hsero_2538; Ref Seq: YP_003775945.1) The H. seropedicae Fnrl, 
Fnr2 and Fnr3 proteins share 38.4%, 37.5% and 26.9% identity 
respectively with E. coli Fnr. As shown in Figure 1, the Fnrl and 
Fnr3 proteins are more similar to each other than to Fnr2. When 
compared with E. coli Fnr, the three Fnr paralogs have characteristic 
sequence features that are hallmarks of Fnr proteins, including three 
conserved N-terminal cysteines plus a central cysteine that are 
thought to co-ordinate the oxygen-labile [4Fe-4S] 2+ cluster 5 . In 
addition all three deduced proteins contain the predicted dimeri- 
zation helix sequence located at the beginning of the C-terminal 
domain, and a helix-turn-helix DNA binding motif characteristic 
of members of the CRP-FNR family. 

A phylogenetic affiliation for the CRP-FNR superfamily of tran- 
scriptional regulators has been proposed 514 . Proteins from the Fnr 
and FnrN groups sense oxygen directly via an [4Fe-4S] 2+ cluster, 
whereas the FixK members, which lack the conserved cysteine resi- 
dues required for ligating the cluster, either sense oxygen indirectly, 
or simply relay the 0 2 signal, via the two-component regulatory 
system FixL-FixJ 13 . In order to classify the H. seropedicae Fnr 



Fnr2_Hs MSKSDNDLSKRLFPVLT CrnCiVKKhCfPHDLSPEET 37 

Fnr_Ec MIP EKRIIRRIQ SGGCAIHCQDCsiSQLClPFTLNEHEL 39 

Fnrl_Hs MNQSCHSTPITTPAN7A7AAPLTHCSSLTAAREAAASSALRSCTACGMHQLCLPMGLDESDM 6 0 

Fnr3_Hs MSSLPASTHKHSPPPVN LHALRASCsaCsMHQLClPMGLDQGDM 44 

: **::*:**.: 

Fnr2_Hs AQFEQLWRRRRVLRGEFAYRAHDECSKIFIVRLGSFKTVRVSRYGGMDVIAFHHTGDLL 97 

Fnr_Ec DQLDNI IERKKPIQKGQTLFKAGDELKSLYAIRSGTIKSYTITEQGDEQITGFHLAGDLV 99 

Fnrl_Hs KRLDKI IGRRK-VARDDFLYRIGDRFTALYAVRVGHFKTYQENLDGDRQITGFQMPGELL 119 

Fnr3 Hs QRLEQVINRRRKVKRDETLYRLNDKFDMLYAIRLGHFKTFQHN PNGGQQITGFQMAGELL 104 



Fnr2_Hs GIEGASHSAYEVDTIALEDSQlCELSFVGLEELSRKIPRLHQQVWRRLSNEVTLMQQQSl 15 7 

Fnr_Ec GFDAIGSGHHPSFAQALETSMVCEIPFETLDDLSGKMPNLRQQMMRLMSGEIKGDQDMIL 15 9 

Fnrl_Hs GMDAISTEQHQCDAVALQDSEVCEIPFARLEQLFGQIPHLLRHFHRIMSHEITSEQNVIM 17 9 

Fnr3_Hs GMDAIGAGHHLCEAVALEDSEVCEIPFASLEDLFRDMPTLLRQFHRMMSLEISREQRVML 



164 



Fnr2_Hs LL-KARSEQRFAFFLLDLSRISSRCGNASTEFSLRMSRTDIGLFLGLTNESMSRLISKFR 216 

Fnr_Ec LLSKKNAEERLAAFIYNLSRRFAQRGFSPREFRLTMTRGDIGNYLGLTVETISRLLGRFQ 219 

Fnrl_Hs LLGNMRAEQRFAAFLVNLSSRYAARGYSSTRFQLRMTRQDVGNYLGLTIES ISRLISKFR 23 9 

Fnr3 Hs TLGSMTAQQKMAAFLLNLSSRYMSRGYSSTRFQLRMTREEIGNYLGLAVESVSRLLTNFK 22 4 



Fnr2_Hs KAGLIDVSVRKVRVLSPSTLEELASGATSWEQLEQQEMNLPCQRNLDLDKRIFNCG 27 2 

Fnr_Ec KSGMLAVKGKYITIENNDALAQLAGHTR NVA 250 

Fnrl_Hs KQ GL LAVE QRDVE VVDLAAL KRLAAGVD ACTATPTRSTS 27 8 

Fnr3 Hs KSGVIEVNHRDVELCDLPTLRAVALGND PCA 255 



Figure 1 | Alignment between H. seropedicae Fnrl, Fnr2 and Fnr3 proteins and E. coli Fnr. Identical amino acids are indicated by asterisks ( * ), 
high similarity amino acids are indicated by colons ( : ) and low similarity amino acids by dots ( . ) . Conserved cysteines required for binding of the 
[4Fe-4S] + are shown in bold and indicated by thick arrows. The double underlined sequence represents the region of the N-terminal sensory domain that 
comprises the eight-stranded (3-roll. The cc-helix required for dimerization is boxed. Highlighted in light-grey is the DNA-binding domain with residues 
that are important for Fnr-box recognition indicated by thin arrows. 
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proteins, we carried out a phylogenetic reconstruction and observed 
that all three are in a clade together with other members of the Fnr 
group of CRP-FNR superfamily from Betaproteobacteria (Supple- 
mentary Fig. SI). This finding is in agreement with the cysteine motif 
arrangement found in the H. seropedicae Fnr proteins, which is char- 
acteristic of the Fnr group and clearly divergent from the FnrN and 
FixK groups commonly represented in the Alphaproteobacteria 5 . 
Among the Betaproteobacteria, H. seropedicae Fnrl, Fnr2 and 
Fnr3 branched into a group with the Fnr proteins from Janthino- 
bacterium sp. Marseille (Minibacterium massiliensis) and H. arseni- 
coxydans. Within this group H. seropedicae Fnr2 is more divergent 
from Fnrl and Fnr3 and also from the /. sp. Marseille and H. arseni- 
coxydans Fnr proteins. 

Construction of H. seropedicae A/nr mutant strains. To study the 
role of the three Fnr proteins, single, double and triple fnr deletion 
mutant derivatives of H. seropedicae strain SmRl were constructed 
using a sacRB::Km cartridge as described in the Methods section. 
This strategy allowed the construction of 7 unmarked deletion 
strains in which all possible combinations of fnr genetic back- 
grounds are available (Supplementary Fig. S2). The promoters and 
non-coding regions of fnr were retained in these ORF deletions, 
enabling transcriptomic analysis of fnr mutant strains. Deletion 
mutants were given the prefix MB, with a number indicating 
which fnr deletion is present (for example MB1 lacks fnrl, whereas 
MB13 lacks fnrl and fnr3 respectively). 

Influence of fnr on growth at low oxygen concentrations. As Fnr 

proteins are known to sense oxygen and have an important role as 
transcriptional regulators during the switch from aerobic to oxygen- 
limiting conditions, we were interested to determine if deletion of the 
three H. seropedicae fnr genes would influence growth under hypoxic 
conditions. Accordingly, we compared the growth curves of wild- 
type strain SmRl with that of the triple deletion fnr strain MB231, 
when grown in malate minimal medium with an initial oxygen 
concentration in the gas phase of 5% and supplemented with 2 or 
20 mM ammonium chloride. Oxygen consumption during growth 
of the cultures was monitored using a gas chromatograph equipped 
with a molecular sieve column and a TCD detector. Under these 
conditions, the rate of oxygen depletion in the gas phase was simi- 
lar for both the wild-type and MB231 strains (Figure 2). However, 
after 4 hours when the oxygen concentration had decreased to 
approximately half of the initial concentration, the growth rate of 
the triple fnr deletion strain was clearly slower than that of the wild- 
type. Moreover, the optical density (O.D 600 ) of the triple mutant 
strain reached only 0.5 and 0.6 after 10 hours growth in the pre- 
sence of 2 and 20 mM NH 4 C1, respectively, compared with an O.D 60 o 
of 0.8 reached by the wild type strain in both ammonium chloride 
concentrations (Supplementary Fig. S3). In contrast we observed no 
difference in growth rate, when both strains were grown under 
aerobic conditions with 20.8% oxygen in the gas phase (Supple- 
mentary Fig. S4). These results imply that the absence of fnr 
imposes a growth rate penalty under oxygen-limiting conditions, 
which may suggest the involvement of at least one of the three Fnr 
proteins in the regulation of terminal oxidases in response to oxygen, 
as observed in other bacteria. 

Transcriptional profiling of wild-type and fnr strains using RNA- 
seq. The influence of the Fnr proteins on global gene expression 
under microaerobic conditions was assessed using RNA sequenc- 
ing. To avoid problems associated with growth rate differences, we 
grew the wild-type and the fnr ablated strain, MB231 under aerobic 
conditions (20.8% oxygen) to an optical density of 0.4 and then 
switched the cultures to microaerobic conditions (initial oxygen 
concentration of 2%) for 1.5 hours prior to RNA extraction. 
Comparison of global gene expression patterns revealed that 187 
H. seropedicae genes were differentially expressed by more than 
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Figure 2 | Influence of fnr genes on growth during oxygen limitation. 

The growth of H. seropedicae SmRl (black squares) and MB231 (grey 
squares) strains were assayed in NFbHP-Malate minimal media 
supplemented with high ammonium concentration (20 mM NH 4 C1) 
under 5% initial oxygen concentration. The oxygen depletion in the gas 
phase was monitored for SmRl (black circles) and MB231 (grey circles). 
Every two hours 0.5 mL samples from the flask gas phase were analysed by 
gas chromatography. The data represents the mean of three independent 
assays performed in duplicate. Error bars indicate standard deviations. In 
some case these are not visible as they are smaller than the graph points. 

3-fold, with p values <0.05. In some cases, depending on the 
genomic context, genes with p-values slightly higher than 0.05 
were also considered as being differentially expressed. Of these, 
143 were down-regulated in the/«r triple mutant strain, indicating 
that these genes are activated either directly or indirectly by Fnr 
under oxygen-limiting conditions. 44 genes were up-regulated in 
the fnr ablated strain, implying that they are targets for Fnr- 
mediated repression. A complete listing of differentially expressed 
genes is provided in Supplementary Dataset 1. 

Of the 187 regulated genes, 70 (37.43%) are classified in the cellular 
process category and 58 genes (31.02%) belong to the metabolism 
category according to the Clusters of Orthologous Genes (COG) 
functional classification 15 (Supplementary Fig. S5a). In the cellular 
process category, 29 genes are related to signal transduction mechan- 
isms, whereas in the metabolism category, 30 genes are related to 
energy production and conversion (Supplementary Fig. S5b). Most 
of the genes from the energy and production subcategory, encode 
important proteins required for synthesis and activity of many of the 
respiratory electron transport chain components (Table 1 and 
Supplementary Dataset 1). These findings suggest that the H. sero- 
pedicae Fnr proteins may facilitate efficient adaptation to the variable 
oxygen concentrations found in different environments. 

H. seropedicae possesses a branched aerobic respiratory chain 
comprising four different types of terminal oxidases 4 . These are 
represented by the aa^-type (cox) and cbb 3 -type (fix) oxidases in 
the heme copper oxidase branch and the bd-type and bo 3 -type oxi- 
dases representing the ubiquinol oxidase branch (Figure 3a). 
Amongst the genes that are potentially activated by Fnr, large 
changes in transcript abundance were observed in genes required 
for the biosynthesis and activity of the cbb } -type heme-copper oxi- 
dase (Table 1). These include thefixNOP operon encoding the struc- 
tural components of this oxidase and the maturation genes fixG, 
fixH, fixl and fixS, which in other bacteria is often organized as an 
operon 16,17 . In contrast, in H. seropedicae the fixG, fixH, fixl and//xS 
genes are dispersed in two distinct operons (Figure 4a), in which fixG 
and Hsero_3199 (fixH) appear to form an operon with the conserved 
putative transmembrane protein Hsero_3198 and fixl and fixS 
apparently form an operon with the heme biosynthesis gene 
hemN 1 " and Hsero_3206, which encodes a conserved hypothetical 
protein (Figure 4a). In agreement with the down-regulation observed 
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Figure 3 | Fnr regulation of components of the electron transport chain (ETC) in H. seropedicae as determined by transcript profiling, (a) Schematic 
representation of the probable organization of ETC branches in H. seropedicae based on the genome annotation, (b) Influence of Fnr on differential 
expression of genes represented in (a). FC indicates fold change. 



in the absence oifnr, both the fixNOP and hemN-Hsero_3206-fixIS 
putative operons have well conserved Fnr-boxes located at positions 
— 143 and — 128 upstream of their respective translational start sites 
(Table 1). The Fnr-Box in the hemN promoter perfectly matches the 
consensus TTGAT-N 4 -ATCAA, while the fixN promoter Fnr-box, 
TTGAT-N 4 -GTCAA, has only one mismatch (underlined) (Table 1). 

Apart from the cbb 3 -type heme-copper oxidase, genes encoding 
the other terminal respiratory oxidases in the H. seropedicae genome 
were not apparently differentially expressed in response to the pres- 
ence of Fnr. Although the genome contains two copies of the coxBA 
operon encoding aa 3 -type oxidases, one of these (Locus Tags: 
Hsero_2311-Hsero_2312) did not appear to be expressed under 
our experimental conditions. The second coxBA operon (Hsero_ 
4160-Hsero_4161) and its associated coxC (Hsero_4157) and coxG 
(Hsero_4159) genes do not appear to be Fnr- regulated. This was also 
the case for the cydAB genes encoding the bd-type oxidase and a bo 3 - 
type oxidase encoded by the cyoABCD operon. However, analysis of 
transcript abundance suggests that expression of several other com- 
ponents of the respiratory chain are subject to regulation by Fnr. 
Transcripts mapping to the petABC operon, which encodes ubiqui- 
nol-cytochrome c reductase (also known as the cytochrome bci com- 
plex or complex III) were down-regulated 7-11 fold in the Fnr 
mutant compared with the wild-type control (Table 1 and 
Figure 3b). In addition, significant differential expression was 
observed for genes encoding cytochrome C551/C552 (Hsero_1104) 
and cytochrome c 553 (Hsero_0153) a c 4 type cytochrome, that poten- 
tially acts as an electron donor to the cbb 3 -type heme-copper oxi- 
dase 19 (Table 1 and Figure 3b). In addition to the apparent influence 
of fnr on the expression of c-type cytochromes (Table 1 and Figure 3), 
the Fnr protein(s) also appear to activate genes required for 



cytochrome c biogenesis. The H. seropedicae genome encodes system 
ITlike machinery for cytochrome c maturation (Supplementary Fig. 
S6), including homologs of the ResB (CcsB) and ResC (CcsA) pro- 
teins, which are proposed to function in the handling of heme and its 
ligation to apocytochrome, and the DsbD (CcdA) and ResA (CcsX) 
proteins that function in the reduction of the disulphide bond in the 
CXXCH heme binding site 20,21 . Transcript profiling revealed that 
genes encoding all four components were significantly down-regu- 
lated in the triple fnr deletion strain (Table 1). The dsbD gene is 
apparently located in an operon downstream of cut A, which is also 
implicated in cytochrome c biogenesis 22,23 . Transcripts mapping to 
the ndh gene, which encodes the non-coupling NADH dehydrogen- 
ase II enzyme, increased in the triple/rcr mutant (Table 1), indicating 
that expression of this enzyme is repressed by Fnr as observed in 
other bacteria 24,25 . Perhaps to compensate for this, the expression of 
Hsero_4284, encoding the energy-conserving NADH dehydrogen- 
ase I is apparently activated 3-fold by Fnr (Table 1). Overall, in 
comparing the global expression pattern of genes involved in elec- 
tron transport, it would appear that Fnr activates expression of genes 
involved in the heme-copper oxidase branch of the aerobic respir- 
atory chain (Table 1 and Figure 3). This involves not only regulation 
of the expression of the cfofc 3 -type oxidase, but also the cytochrome 
bci complex and maturation and expression of c-type cytochromes. 
Hence, Fnr is likely to influence the flow of electrons through the 
cytochrome c branch of the pathway in order to optimize energy 
generation. Additionally, the fnr genes may also control the com- 
position of the quinone pool. Transcripts mapping to ubiF gene 
which encodes a 2-polyprenyl-3-methyl-6-methoxy-l,4-benzoqui- 
nol hydroxylase protein involved in the penultimate step of the ubi- 
quinone biosynthesis pathway and a gene encoding an alternative 
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Figure 4 | Effect of rhri, rhr2 and fer3 mutations on expression of the fixNOP and hemN-Hsero_3206-llxIS operons in H. seropedicae. (a) Schematic 
representation of genomic region encoding the fixNOP and hemN-Hsero_3206-fixIS operons. The genes and predicted functions in the locus are: 
Hsero_3198, transmembrane protein; Hsero_3199, FixH domain containing protein; fixG, iron-sulfur 4Fe-4S ferredoxin transmembrane protein; fixP, 
cbb 3 -type cytochrome c oxidase-subunit lll;fixO, cbb 3 -type cytochrome oxidase-subunit II; fixN,cbb 3 - type cytochrome c oxidase, subunit l;fixS, nitrogen 
fixation protein P-type ATPase protein; fixl, cation transport P-type ATPase protein; Hsero_3206, conserved hypothetical protein; hemN, oxygen- 
independent coproporphyrinogen III oxidase. Black rectangles represent putative FNR-boxes. Genes are not drawn to scale, (b) P-Galactosidase activities 
of fixN::lacZ and hemN::lacZ fusions incubated for 3 hours under the oxygen concentrations of 2.0% (black bars), 4.0% (dark grey bars), 6.0% (light grey 
bars) and 20.8% (white bars). CTRL indicates the SmRl strain carrying the vector plasmid pPW452 (which contains the lacZgene without a promoter). 



2-polyprenyl-6-methoxyphenol hydroxylase (Hsero_4190) also 
required for ubiquinone biosynthesis were significantly down regu- 
lated in the triple fnr deletion strain (Table 1 and Figure 3). In 
common with many other Proteobacteria, genes encoding a respir- 
atory nitrate reductase, organized as a narGHJI-moaA operon, are 
also apparently up-regulated by Fnr in H. seropedicae, as transcripts 
mapping to this operon decreased 16-100 fold in the triple fnr dele- 
tion mutant (Table 1 and Supplementary Dataset 1). In addition, the 
nitrate/nitrate transporter encoded by narKjU operon and the 
nitrate-sensing two component regulatory system narXL are also 
strongly down-regulated in the fnr triple deletion (Table 1 and 
Supplementary Dataset 1). Although this may suggest the potential 
to utilise nitrate as a terminal electron acceptor, the function of this 
respiratory nitrate reductase in H. seropedicae is somewhat enig- 
matic, as various investigators have failed to demonstrate anoxic 
growth of this organism in the presence of nitrate 1,4 . 



Amongst the global changes in transcript abundance, we observed 
differential expression of two of the three H. seropedicae fnr genes 
themselves. Whereas, the transcript abundance upstream oifnr3 did 
not significantly change in the triple deletion mutant, fnrl and fnr2 
were down-regulated 50-fold and 25-fold respectively (Table 1 and 
Supplementary Dataset 1). This suggests several possibilities that are 
not mutually exclusive: (a) Fnr3 is required to activate expression of 
fnrl and fnr2, (b) Fnrl and Fnr2 auto activate their respective pro- 
moters or (c) a combination of Fnr proteins is required to activate 
these promoters. 

Analysis of Fnr regulation of genes encoding the cbb 3 -type respi- 
ratory oxidase. To confirm the involvement of Fnr in co-regulation 
of the fixNOP and hemN-Hsero_3206-fixIS operons we constructed 
fixN::lacZ and hemN::lacZ transcriptional fusions and analysed the 
expression of P-galactosidase in the various fnr mutant strains 
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Figure 5 | Fnr proteins influence the cytochrome content in H. seropedicae. (a) Bacterial suspensions of H. seropedicae SmRl and fnr mutant strains. 
Cells, from 100 mL cultures, were collected by centrifugation and resuspended in 10 mL of buffer (100 mM NaCl and 50 mM Tris.HCl pH 7.5). 
(b) Reduced minus oxidized visible absorption spectra of protein extracts from H. seropedicae SmRl and fnr mutant strains. For simplification, only the 
data for SmRl (blue), MB1 (red), MB2 (green), MB3 (black) and MB231 (orange) are shown, (c) Heme stained gel of H. seropedicae protein extracts. 
Samples (50 ug protein per lane) from H. seropedicae SmRl and fnr mutant strains were separated by 12% Tris-Tricine SDS-PAGE and stained for 
covalently bound heme with o-dianisidine. On the left the heme stained bands are labeled as 1, 2, 3, 4 and 5 from the top to the bottom of the gel. The 
apparent molecular masses of proteins (kDa) are indicated on the right. MW: Spectra™ Multicolor Broad Range Protein Ladder (Fermentas). The strains 
MB13, MB21 and MB23 gave similar spectra and heme stain profiles to MB231, MB1 and MB3, respectively (Supplementary Fig. S8). 



compared with the parental strain. Expression of the fixN::lacZ and 
hemNr.lacZ fusions is apparently regulated by oxygen levels in the 
wild-type background since we observed a reduction in promoter 
activities upon exposure to increasing oxygen concentrations 
(Figure 4b). Consistent with the transcriptomics data (Table 1), 
both operons are apparently subject to regulation by Fnr. Notably, 
expression from the fixN and hemN promoters was significantly 
reduced in strains that lack either fnrl or fnr3, but activity was 
equivalent to the parental strain in the fnrl deletion strain MB2 
(Figure 4b) . This implies that both Fnrl and Fnr3 are required to 
activate expression of the fixNOP and hemN-Hsero_3206-fixIS 
operons and that Fnr2 is not involved in the regulation of 
expression of the cbb 3 -type respiratory oxidase. Accordingly, under 
oxygen -limiting conditions, no growth penalty is observed for the 
MB2 strain (Supplementary Fig. S7). 

Fnr influences the cytochrome content of H. seropedicae. As the 

transcriptome analysis implicates Fnr as a regulator of genes 
involved in cytochrome c biogenesis, we compared the spectral 
features of wild-type and fnr deletion strains. We noticed that 
strains lacking the fnrl gene were deficient in a pink pigment 
when cultured in liquid media (Figure 5a). To further explore this 
observation we analysed reduced minus oxidized spectra of protein 
extracts obtained from SmRl and the fnr deletion strains (Figure 5b 
and Supplementary Fig. S8). Spectra of the wild type strain were 
consistent with the presence of c-type (a-band located around 
550 nm) and fo-type (a-band shoulder around 560 nm) cyto- 
chromes in the protein extract. Similar spectral features were 
found in strains lacking either fnrl or fnr3. However, all strains 
lacking fnrl appeared to be deficient in cytochrome content, which 



may account for the observed differences in culture pigmentation. In 
order to obtain further biochemical support for the spectral features 
observed, we stained protein extracts from SmRl and fnr mutant 
strains for covalently bound heme (Figure 5c and Supplementary 
Fig. S8). In the wild-type strain SmRl, we detected five bands, 
which presumably represent c type cytochromes. The protein of 
approximately 34 KDa (band 1) could represent FixP by 
comparison with the heme staining profile of the fixN mutant 
strain, RAM21 (Supplementary Fig. S9) and by analogy with 
studies on the FixNOQP proteins from other bacteria 26,27 . Notably, 
the level of this protein was significantly diminished in strains 
lacking fnrl, consistent with decreased expression of the fixNOP 
operon observed in the transcriptome (Table 1 and Figure 3) and 
ZacZ-fusion analysis (Figure 4). The bands 2, 4 and 5 may represent 
PetC, cytochrome c 553 (Hsero_0153), and cytochrome c 551 /c 552 
(Hsero_1104), respectively, based on the apparent molecular 
masses and expression pattern of these proteins, which were 
identified as being activated by Fnr in the transcriptome analysis 
(Table 1 and Figure 3). All five c-type cytochromes, including band 
4, were absent in strains lacking both fnrl and fnr 3 (MB 13 and 
MB231). This is in agreement with the loss of the cytochrome ot- 
band in the UV-visible difference spectra in strains lacking both the 
fnrl and fnr3 genes (Figure 5b and Supplementary Fig. S8). Taken 
together with the data from transcription profiling, these results 
suggest that both Fnrl and Fnr3 are necessary to maintain the 
level of c type cytochromes in H. seropedicae under microaerobic 
conditions. 

Deletion of the three fnr genes impairs nitrogenase activity and 
growth on dinitrogen. In the analysis reported so far, strains were 
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Figure 6 | Influence of Fnr and FixN on nitrogenase activity. The acetylene reduction assay was performed as described in Methods using strains grown 
in (a) semi-solid medium supplemented with 0.5 mM of sodium glutamate or (b) liquid medium supplemented with 0.5 mM of ammonium chloride 
under 5% initial oxygen in the gas phase. In (b) samples were taken from the culture to measure the growth curve (primary y axis) of SmRl (black 
squares), RAM 21 (dark grey triangles) and MB231 (light grey circles). Black, dark grey and light grey bars (secondary y axis) indicate the nitrogenase 
activity of SmRl, RAM21 and MB231, respectively. Data represent the average of two independent experiments performed in duplicate. Error bars 
indicate standard deviations. 



grown in minimal media containing a high concentration of fixed 
nitrogen, which represses nitrogen fixation in H. seropedicae. Since 
the cbb^-type heme-copper oxidase is known to have an important 
role as a terminal oxidase that supports nitrogen fixation under 
microaerobic conditions in the Rhizobacteriaceae 2S ~ io and is subject 
to regulation by Fnr proteins, we were interested to determine if 
nitrogen fixation is influenced by the presence of the Fnr paralogs 
in H. seropedicae. 

When fnr mutant strains were grown under N-deficient condi- 
tions (with 0.5 mM sodium glutamate) in semi-solid medium and 
tested for the ability to reduce acetylene as a measure of nitrogenase 
activity, no significant differences were observed in comparison with 
the wild-type strain (Figure 6a). The RAM21 strain (fixN mutant) 
was also not deficient in acetylene reduction when grown under these 
conditions (Supplementary Fig. S10). However, since semi-solid 
medium enables bacteria to move towards optimal oxygen concen- 
trations appropriate for growth, we sought a more rigorous method 
to determine the influence of limiting oxygen on nitrogenase activity 
in the mutant strains. When the fixN mutant strain RAM21 was 
grown in nitrogen-deficient liquid medium, under conditions of 
oxygen limitation (initial oxygen concentration of 5% in the gas 
phase) we observed that acetylene reduction was not severely com- 
promised in comparison with the wild-type (Figure 6b). This sug- 
gests that the cfr£> 3 -type oxidase is not required to support 



nitrogenase activity in H. seropedicae under the oxygen-limiting con- 
ditions imposed in this experiment. In contrast, the nitrogenase 
activity of the MB231 mutant strain was severely impaired compared 
to the wild-type strain under these conditions (Figure 6b). Given the 
pleiotropic effects on expression of the electron transport compo- 
nents in the triple fnr deletion strain it is likely that the electron flux is 
insufficient to support nitrogenase activity in the fnr deletion strain 
under oxygen-limiting conditions. In agreement with this result, we 
observed that diazotrophic growth of the triple/nr mutant strain was 
also compromised in liquid N-free medium (Supplementary 
Fig. SI 1). 

Discussion 

In order to survive in rapidly changing environments and explore 
diverse habitats, many bacteria adjust the composition of their res- 
piratory chains to cope with fluctuating oxygen concentrations. In 
many cases this involves regulation of the expression of terminal 
oxidases in order to optimize energy generation. The global tran- 
scriptional regulator Fnr, and its various orthologs, provide a wide- 
spread mechanism for sensing oxygen and communicating this to 
the transcriptional apparatus in order to balance the levels of differ- 
ent terminal oxidases, according to prevailing environmental condi- 
tions. Accordingly, we observe that Fnr is required to activate the 
expression of genes required for the synthesis and activity of the 
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high-affinity cbb 3 -type heme copper oxidase in H. seropedicae, as is 
the case in other Proteobacteria 16,31 . To balance respiratory require- 
ments under oxygen-limiting conditions, Fnr and its orthologs com- 
monly participate in negative regulation of the expression of other 
terminal oxidases, for example, the bd-type and bo 3 -type oxidases in 
E. coli and A. vinelandii 32,33 and the bo 3 -type and CIO oxidases in 
Pseudomonas putida 34 . However, our transcriptomics data indicate 
that this is not the case in H. seropedicae. We only observe differential 
expression of the genes encoding the cbb 3 -type oxidase. Expression 
of the other terminal oxidases is not apparently affected by absence of 
the three H. seropedicae Fnr proteins. However, in contrast to other 
well-studied systems, the Fnr proteins in H. seropedicae appear to 
have a major influence on the composition of the complete electron 
transport chain that feeds electrons from NADH, through the ubi- 
quinone pool to the cytochrome bcj complex and onto the c-type 
cytochromes that are substrates for the heme-copper oxidases. This 
is clearly demonstrated by the depletion of c-type cytochromes and 
the down regulation of genes encoding the various components of 
this branch of the electron transport chain in the triple fnr mutant. 
Hence it would appear that the Fnr proteins play a major role in 
regulating the configuration of the H. seropedicae electron transport 
chain in order to exploit respiratory flexibility and optimize energy 
coupling in response to oxygen availability. 

The cbb 3 -type heme copper oxidase is likely to be required for 
growth at very low (<0.5%) oxygen concentrations 35 and can sup- 
port symbiotic nitrogen fixation at nanomolar levels of dissolved 
oxygen (reviewed in 36 ). Nevertheless, the fixN insertion mutant of 
H. seropedicae RAM21 was competent to support nitrogenase activ- 
ity under oxygen-limiting conditions implying that the cbb 3 -type 
oxidase is not required to support nitrogen fixation in this organism. 
Perhaps this result is not surprising, given that, to our knowledge, 
this oxidase is not required for nitrogen fixation in other free-living 
diazotrophs. By analogy with other nitrogen-fixing bacteria, it is 
possible that the bd-type oxidase supports nitrogenase activity in 
H. seropedicae. This oxidase is critical for microaerobic diazotrophy 
in Klebsiella pneumoniae 37 , it provides respiratory protection for 
nitrogenase in Azotobacter vinelandii 3 "' 39 and it is utilized as a ter- 
minal oxidase to support symbiotic nitrogen fixation in Azorhizo- 
bium caulinodans 35 . In contrast, the triple fnr deletion strain was 
compromised with respect to both nitrogenase activity and diazo- 
trophic growth, which presumably reflects the major role played by 
Fnr in reconfiguring the electron transport chain under oxygen lim- 
iting conditions in H. seropedicae. 

The results presented here do not provide a rationale for the exist- 
ence of multiple Fnr proteins in H. seropedicae. Potentially, each 
ortholog may exhibit differential sensitivity to oxygen, recognize 
different DNA targets or have different propensities to dimerise 
under aerobic conditions. The transcript profiling reveals that 
expression oifnri is constitutive, whereas fnrl and fnr2 expression 
is apparently positively controlled by one or more of the Fnr ortho- 
logs. The. fnrl gene, which is located close to the genes required for 
the maturation and activity of the cbb 3 -type oxidase, appears to play a 
critical role in regulating transcription of the fixNOP and hemN- 
Hsero_3206-fixIS operons and in controlling the expression of c-type 
cytochromes. Although, fnr3 also appears to be required to express 
the cfofr 3 -oxidase, we cannot rule out the possibility that it is required 
to activate transcription of fnrl. In contrast, fnrl does not appear to 
be required, either for expression of this oxidase or c-type cyto- 
chromes, under the experimental conditions employed here. 
Further detailed characterization of the three Fnr paralogs will be 
necessary in order distinguish their precise roles in gene regulation in 
H. seropedicae. 

Methods 

Bacterial strains and plasmids. H. seropedicae and E. coli strains and plasmids used 
are listed in Supplementary Table Si. 



Growth conditions. E. coli strains were grown at 37°C in LB medium 40 . H. 
seropedicae strains were grown at 30°C in NFbHP-Malate medium 41 supplemented 
with NH4CI or 0.5 mM sodium glutamate. Appropriate antibiotics were used when 
required. For experiments requiring different oxygen concentrations, the air in the 
gas phase of Suba Seal® stoppered culture flasks was exchanged by injecting argon 
into the flasks for 30 minutes. To obtain different oxygen levels a given volume of air 
was injected back into the flask. The oxygen levels in the gas phase were verified by gas 
chromatography using a molecular sieve column and a TCD detector. 

Identification and Analysis of H. seropedicae Fnr Orthologs. The sequences of 
Fnrl, Fnr2 and Fnr3 from H. seropedicae SmRl were aligned with E. coli K12 substr. 
MG1655 Fnr (Ref seq: NP_415850.1) using Muscle software 42 . The presence of 
conserved domains in the Fnr homologs was investigated by submitting sequence 
search to the Pfam database {http://pfam.sanger.ac.uk/) 43 . 

Phylogenetic Analysis. Amino acid sequence retrieval was performed by using a 
BLASTP search 44 against the nonredundant NCBI database. The proteins used as 
queries were Fnrl, Fnr2 and Fnr3 from H. seropedicae, (Ref seq: YP_003776587.1, 
YP_003775788.1 and YP_003775945.1, respectively) FixKl and FixK2 from B. 
japonicum (Ref seq: NP_772701.1 and NP_769397.1, respectively), Fnr from E.coli 
K12 substr. MG1655 (Ref seq: NP_415850.1) and CydR from A. vinelandii (Ref seq: 
YP_002799173.1). A limited number of sequences were selected, eliminating 
redundant information, while maintaining representative taxonomic diversity. All 
sequences selected (Supplementary Table S2) were checked to have the CRP-FNR 
superfamily protein signatures, CNMP_Binding _3 PS50042 and HTH_CRP_2 
PS51063, proposed by PROSITE 45 . For phylogenetic tree reconstruction, an amino 
acid alignment was made using Muscle 42 . Maximum likelihood (ML) trees were 
derived using the JTT matrix-based model 46 after bootstrapping 1,000 replicates of 
each original data set 47 using the MEGA 5.05 software 48 . 

Construction of H. seropedicae SmRl fnr deletion and fixN insertional mutant 
strains. An allelic exchange strategy was used to generate derivatives ofH. seropedicae 
SmRl containing/hr orthologs deletions (Supplementary Fig. S2) and a tetracycline 
resistance cassette insertion intofixN (Supplementary Fig. S12). The primers used in 
this work are show in Supplementary Table S3. For construction of the allele exchange 
plasmids for fnr mutation the upstream and downstream regions of fnrl, fnrl and 
fnr3 were amplified by PCR. These fragments were then ligated to generate fnr 
deletions, which were cloned into Hindlll and BamHI sites of the suicide plasmid 
pSUP202, to generate pMBBID, pMBB2D and pMBB3D (Supplementary Table SI). 
The nptl-sacB-sacR cartridge (from pMH1701) was then inserted into BamHI site. 
We generated three suicide plasmids: pMBBlDS for the 279 bp deletion of fnrl, 
pMBB2DS for the 276 bp deletion of fnrl and pMBB3DS for the 267 bp deletion of 
fnr3 (Supplementary Table Si). 

Conjugation was performed between E. coli S17.1 containing the plasmid of 
interest and H. seropedicae recipient strains. Conjugation was performed on NFbHP- 
Malate/LA (3 : 1) agar by mixing recipient and donor strains in two proportions (50/1 
and 10/1). Transconjugants were selected on NFbHP-Malate agar supplemented with 
20 mM NH4CI and antibiotics. One mutant strain resulting from single cross-over 
was grown overnight in liquid NFbHP-Malate plus 20 mM NH 4 C1 without anti- 
biotics at 30°C, After incubation 250 uL of the culture were plated on NFbHP-Malate 
agar, supplemented with 20 mM NH 4 C1, 5% sucrose, 5 ug/mL nalidixic acid and 
80 jig/mL streptomycin. Sucrose is toxic to bacteria that express the sacB gene, 
therefore only strains that lost the sacRB-Km R cassette by a second homologous 
recombination event could grow under these conditions. The mutant strains were 
analysed by PCR using primers (Supplementary Table S3) external to fnrl, fnrl and 
fnr3. To construct double and triple H. seropedicae mutant strains the process 
described above was repeated using different allele exchange plasmids. 

For fixN mutagenesis, the plasmid pHSl7058Hll (Genopar Consortium) 4 con- 
taining the/ixJVgene was subject to a transposition reaction using the EZ-Tn5 <TET- 
1 > Insertion Kit (Epicentre) . After confirmation of the TET- 1 insertion into the fixN 
gene (pHS17058HllTc), a cloramphenicol cassette from pTnMod-OCm plasposon 
was inserted into the plasmid, outside thefixN gene, to yield the plasmid 
pHSl7058HllTcCm, which was then electro -transformed into H. seropedicae SmRl. 
A transconjugant resistant to tetracycline and sensitive to cloramphenicol was 
selected and named RAM2 1 . The double recombinant was confirmed by DNA 
hybridization (Supplementary Fig. S12). 

RNA isolation and RNAseq library construction. For total RN A extraction, we grew 
H. seropedicae SmRl (wild-type) and MB231 strains (triple/«r mutant) under aerobic 
conditions to an optical density of 0.4 (cultures were shaken at 120 rpm in air) and 
then the cultures were switched to microaerobic conditions (initial oxygen 
concentration 2%) for 1.5 hours. After collection of the cells by centrifugation, the 
total RNA was isolated using the RNA RiboPure 1M -Bacteria kit (Ambion). Two 
rounds of enrichment were performed using MICROB£xpres5 1M kit (Ambion) for 
removing ribosomal RNA from total RNA samples. Approximately 200 to 500 ng of 
mRNA enriched RNA were used for cDNA synthesis and library construction using 
the SOLiD 1M Whole Transcriptome Analysis Kit (Life Technologies). The libraries 
were amplified by emulsion PCR using SOLiD ePCR kit and sequenced in a SOLiD 4 
System (Life Technologies). 

Read Mapping, Differential Expression Analysis and Fnr Binding Site Prediction. 

The reads were mapped against the H. seropedicae SmRl genome (NC_014323) as 
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reference using the CLC Genomics Workbench package. Read counts table was 
exported into the RobiNA software 49 and both normalization and statistical 
evaluation of differential gene expression was performed by DESeq 50 with ap-value 
cut-off of 0.05 using the Benjamini-Hochberg method for multiple testing 
correction 51 . Genes with fold change lower than three were excluded from the 
analysis. Genes with fold change lower than 3 orp-value slightly higher than 0.05 were 
included in the analysis when genome context indicated they are part of operon with 
genes differentially expressed according to the above parameters. 

Potential Fnr-binding sites were located using the sequence motif search facilities 
of PEPPER 52 (http://pepper.molgenrug.nl/index.php) and Virtual Footprinting 53 
(http://prodoric.tu-bs.de/vfp/vfp_promoter.php) using default search parameters. 

Construction of transcriptional fusions. To study the expression of fixNOP and 
hemN-Hsero_3206-fixIS operons, the putative promoter regions were amplified by 
PCR, cloned into pT257R/T and then subcloned into PstI and Bglll sites of pPW452 
to generate the plasmids pPWPFN (fixN::lacZ fusion) and pPWPHN (hemNr.lacZ 
fusion) (Supplementary Table Si). All constructs were verified by DNA sequencing. 

p-Galactosidase activity. To analyse the activities of fixNr.lacZ and hemNr.lacZ 
fusions, strains were grown to an O.D 600 = 0.3, collected and then resuspended in 
NFbHP Malate liquid medium supplemented with 20 mM NH 4 C1 to an O.D 600 = 
0.05. The cells were incubated for 3 hours either under air or at initial concentrations 
of 2, 4 or 6% oxygen. After incubation samples were taken for |3-galactosidase activity 
determination as described 54 . The results are expressed in Miller units (MU). 

Acetylene reduction assay. The acetylene- reduction assay was used to determine 
nitrogenase activity on free-living cultures 55,56 . Ethylene formation was determined 
either by using a Varian Star 3400 CX gas chromatograph equipped with a Porapak N 
column. 

Freshly grown cultures were used for inoculating NFbHP -Malate semi solid 
medium (0.17% agar) containing 0.5 mM sodium glutamate followed by eighteen 
hours incubation at 30 C. Acetylene (10%) was injected and after incubation for one 
hour at 30 C, 0.5 raL samples were collected for determination of produced ethylene 
by gas chromatography. The same procedure was used for assaying nitrogenase 
activity in liquid media, except that the cultures were collected by centrifugation 
(3220 g, for 4 min at room temperature), re-suspended in NFbHP-Malate supple- 
mented with 0.5 mM of NH 4 C1 to an O.D of 0.2. The cells were incubated at 5.0% 
initial oxygen concentration in the gas phase as described. Nitrogenase activity is 
reported as nmol of C 2 H 4 produced per minute per mg protein. Whole cell protein 
concentration was determined by the Bradford method 57 after overnight lysis with 

0. 2 mM NaOH. 

Reduced minus Oxidized Spectra and Heme Stain. For preparation of protein 
extracts H. seropedicae strains were cultured in 250 raL erlenmeyer flasks containing 
50 mL of NFbHP-Malate supplemented with 20 mM NH 4 C1. After overnight 
incubation, the cells were collected by centrifugation and re- suspended in 2 mL of 
sonication buffer (100 mM NaCl and 50 mM Tris.HCl pH 7.5). Cells were broken by 
sonication, and cell debris was separated from the protein extract by centrifugation 
(6000 g, 30 min). The supernatant was collected and used for further analysis. 

UV -visible difference spectra of 0.3 mg.mL -1 protein extracts were recorded in 
1 cm path-length quartz cuvettes at room temperature on a Shimadzu UV-2501 PC 
spectrophotometer. Reduced-minus-oxidized difference spectra of the protein frac- 
tions were recorded by measuring the dithionite-reduced spectrum of the sample 
against the air- oxidized one. 

For the covalently bound heme stain, 50 ug of the protein extract, prepared as 
described above, was loaded without boiling onto a 12% Tris-Tricine SDS-PAGE 
gel 58 , and stained using o-dianisidine 59 . 
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